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Depletion of background galaxies due to the cluster lens 
CL0024+1654: U and R band observations 



o: 
o 

O' 
(N; 

c/3 



6. E. Rognvaldsson, 1 * T. R. Greve, 2 J. Hjorth, 2 

E. H. Gudmundsson, 3 V. S. Sigmundsson, 3 P. Jakobsson, 2 ' 3 

A. O. Jaunsen, 4 L. L. Christensen, 5 E. van Kampen 6 and A. N. Taylor 6 

1 NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

2 Astronomical Observatory, Univeristy of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark 

3 Science Institute, University of Iceland, Dunhaga 3, IS- 107 Reykjavik, Iceland 

4 Institute of Theoretical Astrophysics, Postboks 1029 Blindern, N-0135 Oslo, Norway 

5 ST-ECF/ESO, Garching, Germany 

6 Institute for Astronomy, Blackford Hill, Edinburgh EH9 3HJ, UK 



> 

in 

^- 

o 
o\ 
o 
o 
o 

S3 

Oh: 

6 

c3 



X 



Accepted . . . Received 2000 July 17; revised 2000 August 30 



ABSTRACT 

We have obtained U and R band observations of the depletion of background galaxies 
due to the gravitational lensing of the galaxy cluster CL0024+1654 (z = 0.39). The 
radial depiction curves show a significant depletion in both bands within a radius of 
40" — 70" from the cluster center. This is the first time depletion is detected in the 
U band. This gives independent evidence for a break in the slope of the U band lu- 
minosity function at faint magnitudes. The radially averaged R band depletion curve 
is broader and deeper than in the U band. The differences can be attributed to the 
wavelength dependence of the slope of the luminosity function and to the different 
redshift distribution of the objects probed in the two bands. We estimate the Ein- 
stein radius, r&, of a singular isothermal sphere lens model using maximum likelihood 
analysis. Adopting a slope of the number counts of a = 0.2 and using the background 
density found beyond r = 150" we find r E = 17" ± 3" and r E = 25" ± 3" in the U 
and R band, respectively. When combined with the redshift of the single background 
galaxy at z — 1.675 seen as four giant arcs around 30" from the cluster center, these 
values indicate a median redshift in the range (zs) ~ 0.7 to 1.1 for the (7ab> 24 mag 
and i?AB> 24 mag populations. 

Key words: Gravitational lensing - galaxy counts - galaxies: clusters: general - 
galaxies: clusters: individual: CL0024+1654 - ultraviolet: galaxies 



1 INTRODUCTION 

Gravitational lensing by clusters of galaxies affects the ap- 
parent distribution of background galaxies on the sky, pro- 
vided that the slope of the galaxy number counts differs from 
a critical value corresponding to a balance between the di- 
lution due to local stretching of the background sky and the 
density enhancement due to magnification of faint sources 
above the observed magnitude limit. Normally, the slope 
levels off at faint magnitudes, causing an apparent deple- 
tion of galaxies behind strong gravitational lenses relative to 
the field. This depletion signal is, under fairly general condi- 
tions, directly related to the magnification of the lens , which 
in turn may be used to find its surface mass density. Broad- 



D 



hurst et al. (1995) first pointed out the potential of this as- 
pect of the magnification bias to estimate cluster masses and 
density profiles, using the so-called number count method. 
This method was fi rst used for absolute mass estimation by 
Taylor et al. (1998), who also discuss various observational 



issues. Models for relating the magnification to the surface 
mass density (or convergence), going b eyond the weak lens - 
ing approximation, were considered by van Kampen (199c) 
and Dye & Taylor (1998). Such models, which take into ac- 



count the presence of a shear term in the magnification, 
are of great importance in combination with the number 
count method, since the shear patter n can not be measured 
from the background densities alone. Schneider et al. (200C) 
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have compared the number count method to shear measure- 
ments, using maximum likelihood analysis. They find that 
the number count method does a better job of estimating 
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cluster density profiles, provided that the unlens ed number 
density of g alaxies j s known to a good precision. Bezecourt 



stein radius of CL0024+1654, inferred from the position of 
giant arcs around 30" from the cluster center, and pr evi- 



et al. (199S ) have investigated the wavelength dependence of ous de tection of depletion in the I and B band (Fort et al 



depletion curves, while Mayen fc Soucail (2000]) investigate 
how depletion curves depend on the cluster density profile. 
They also investigate how the filter bands used in obser- 
vations affect the results through different sampling of the 



(1997)) make this cluster a good candidate for depletion 



backgrc und population. Gray et al. (200C|) show that the ef- 
fects of incompleteness on mod el parameters estimate d with 
maximum likelihood methods ( Schneider et al. (200C )) may 
be neglected, at least under the standard assumption that 
the intrinsic luminosity function is a power law. 

A few other methods based on gravitatio nal lensing 



studies. Moreover, its mass and mass distribution have been 
studied carefully with other methods. 

F r om the redshift catalo g co mpiled by 
(19850 , pchneider et al. (19860 , and~ 



Dressier . 



a velocity dispersion of a v = 1200 km s 1 
suiting in a mass of roughly 7 x 10 hTg 1 M© 



Dressier et a l 
Gunn (1992^ " 



is found, re- 
within a ra- 



have been applied to cluster mass estimation. Dye et al 



(2000hhave studied the effects of the lens magnification on 
the luminosity function of background population, and ap- 
plied this method to Abell 1689. Weak shear analysis, where 
the average distortion of the shape of background galaxies is 
used to estimate the shear fi eld of the lens, has been a p plied 
to a nu mber of cluste rs (e.g. Smail fc Dickinson (199£ ) ; [Lup- 
pino fe|Kaiser (1997); |ciowe et al. (199g~Hoekstra et al 



dius of 0.5h^ Mpc. (The Hubble parameter is Ho = 50/iso 
km/s/Mpc.) A much larger sample of redshifts in the cl uster 
field has recently bee n obtained by Dressier et al. (199£) and 
Czoske et al. (1999 ). Analysis of this data leads to a con- 
siderably lower velocity dispersion, a ~ 700 km s , for the 
227 cluster objects identifie d. The estimated m ass is conse- 
quently lower. For example |Czoske et al. (1999) find a mass 



(200 0pJo"ffre et al. (2000[ ); |Athreya et al. (20000 , and refer- 
ences therein). Strong lensing effects have also been used to 



probe t he density distribution in clusters (e.g. Tyson et al 
(19981) 



of 1.4 x lO 14 /^ 1 M Q within 0.5/i^j 1 Mpc. 

From a detailed inversion of the lens using mainly the 
giant arcs, whic h are the multiple i mages of a single back- 
ground galaxy, Tyson et al. (199£) find the central mass 
distribution to be very relaxed, with no signs of substruc- 
ture within 220/i^q 1 kpc. On the other hand, Broadhurst 



et al. (200C) find evidence for substructure on the cluster 

l 



Broadhurst et al. (200C)). and statistical studies of core scale. Their mass estimate for the innermost 220/i 



arclets are being pursued (e.g. Bezecourt et al. (1999 )). The 
theoretical and observational aspects of gravitational lensing 
and its app lications have recently been rev iewed in 
(1999Qind |Bartelmann fc Schneider (1999J ). 



kpc is 2.6 x 



io 14 ^ 1 



Mp, , or somewhat 



Mellier 



3.3 x W h~ L M Q estimated by Tyson et al. (199£ 



lower than the 
Their 



To date, depletion of background galaxies h as been ob- 
served in connection with a handful of clusters. Fort et al _ 
(19970 vere the first to report on depletion in the B and I 



model does however predi ct an unrealistically hi gh velocity 
dispersion for the cluster ( ghapiro fc Iliev (200C )). 

On a much larger scale, Bonnet et al. (1994|) using weak 



shear analysis find a mass of 4 x 10 h 50 M@ within a radius 



of 3fo Rfl Mp c. An alysis of X-ray obse rvations by Bohringer 



band b jhind CL0024+1654. Their observations were how- 



ever plagued by a low density of obje cts and an uncertai n 
background estima te as discussed by van Kampen (1998). 
Taylor et al. (1998J ) found depletion in the / band in Abell 
1689, after careful removal of cluster and foreground oh- 



et al. (2000] ) and Soucail et al. (2000) shows a discrepancy 
of about a factor 1 — 3 between X-ray and lensing mass 
similar to what is found in many other clusters. The differ- 
ence is most probably attributable to simplified as sumptions 



about the dynamical state of the cluster. Finally, yan Kam- 



iag; additional color information. [Maven & Soucail| pen (19981) used the depletion observations by 



Fort et al 



jects m 

(2000[)~ ind depletion of faint B, V, R, and J objects towards (1991) to constrain the mass of the cluster, noting that a bet- 



MS1008-1224. This cluster can possibly be used to investi- 
gate how background clustering affects the depletion anal- 
ysis, since Athreya et al. (200CJ ) have tentatively detected 
another cluster i n the background (at z — 0.9). Finally, 
Gray et al. (200C ) have explored depletion at near-infrared 
wavelengths in Abell 2219. Depletion analysis in the infrared 
(IR) benefit from a low slope of the number counts already 
at moderately faint IR magnitudes, but the low background 
number density may outweigh the benefits. 

In our work, we apply the number count method in the 
U and R band to the rich cluster CL0024+1654 (z=039) 



Koo 



first observed to be a strong gravitational lens by 
(19880. [The U filter is well suited for depletion studies, pro 
vided that the observations are deep enough to probe galax- 
ies fainter than Uab~ 25.75 where th e slope of the luminos - 
ity function in the U band levels off (Williams et al. (1996); 



Pozzettiet al. (19980 , but see e.g. [Hogg et al. (1997J ); [Volon, 
teri et 1 1. (200C ) for evidence against this) , while the slope in 
the R band is low enough to accommodate depletion studies 
given sufficient number of background objects. With suffi- 
ciently deep observations, depletion analysis in these two 
widely separated bands will also reflect the different sam- 
pling of the background redshift distribution. The large Ein- 



ter estimate of the unlensed background density was needed 
to give reliable results. 

In order to constrain the mass distribution of the cluster 
using our data, additional information is needed on various 
properties of the background sources. The unlensed average 
density, redshift distribution, and the slope of the number 
counts of the background sources all play a role in relating 
the observed depletion to the convergence of the lens. In this 
paper, we present the first U band observations of depletion 
and apply simple parametric modeling to our U and R band 
data to investigate the dependence of the model on the a 
priori parameters. A more detailed investigation of the mass 
distribution implied by our data is being undertaken (van 
Kampen et al., in preparation). 

The rest of this paper is organized as follows. Basic 
description of the data is given in Section H. The object 
extraction and analysis of the resulting catalogs is discussed 



in Section 
in Section 



and the results are discussed and summarized 



2 THE DATA 
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2.1 Data acquisition and reduction 

The data were obtained with the ALFOSC instrument on 
the Nordic Optical Telescope in August 1999. The field-of- 
view was 6.5D' with a pixel size of 0.188". The total inte- 
gration time in the (Cousins) U band was 37 ksec and 8.7 
ksec in the (Cousins) R band with average seeing of 1.1" and 
1.0", respectively. The formal 3<r detection limits in the final 
drizzled images are UAB,lim— 27.2 mag and -RAB,iim= 27.5 
mag. (We use AB magnitudes throughout, using a simple 
shift of +0.71 mag and +0.199 mag to convert fro m the ob- 
served Cous ins U and R magnitudes, respectively ( Fukugita 
et al. (jgjlj)). 

Standard bias and flat field corrections were made to the 
data before stacking t he individual frames usi n g tools from 
the dit herll p a ckage (Fruchter fc Hook (1997J ); |Mutchlei~& 
Frucht j~(l"997T ); [Hook et al. (1999D ) in IRAF. Shifts were es- 
timated by cross-correlating images and cosmic ra ys masked 
out usi ng the blo tting technique as described by Mutchler 
& Fruc'iter (1997). The original pixel scale was retained and 
no attempt was made to correct for field distortions. 



2.2 Photometry and completeness 

The c alibration aga inst photometric standards (m ainly in 



Landolt 



M92 (Davis (1999)) and a few standards from 
(1983)) gives formal errors of ctab = 0.02 mag and <tae/ = 
0.05 mag. Typical photometric errors estimated by the de- 
tection software are between 0.1 — 0.3 mag for the magnitude 
ranges 24 < (R, U)ab < 27. (This is the magnitude range 
used below to select background objects for the depletion 
analysis.) 

In order to estimate the level of completeness as a func- 
tion of magnitude in our data, 300 artificial objects were 
added to the science frames using the artdata task in IRAF. 
Completeness was estimated as the ratio between detected 
synthetic objects and the number of synthetic objects added. 
In order to improve the statistics this procedure was re- 
peated 100 times, each time adding objects at random posi- 
tions. The final average completion curves indicate that our 
sample is 50% complete at Uab= 26.7 mag and Rab= 26.5. 

The fraction of detections due to spurious noise peaks 
will steadily increase towards faint flux levels. To investi- 
gate this, frames containing only noise wer e generated and 
using SExtractor (Bertin & Arnouts (1996)), 700 apertures 



were dropped onto the frames. To ensure that the noise-only 
frames would mimic the noise characteristics of the science 
frames as closely as possible, they were convolved with the 
noise correlation function found in the science frames. By 
counting the number of detections per magnitude bin we 
can readily quantify the contamination from noise peaks. 
It was found that only 7% of the detections have U mag- 
nitudes brighter than 27 mag while in the R frame about 
9% have magnitudes brighter than 27 mag. With this low 
level of noise contamination, we decided not to make any 
corrections for false detections. 



3 ANALYSIS 

The number count method is based on the following prop- 
erty of cluster lenses. Under the standard assumption that 



the number density of background galaxies follows a power 
law, no(S) ex S'"' 3 for a given flux limit S, one finds a sim- 
ple relation between the magnification of the lens, n, and 
n(6; > S), the number density of background objects as ob- 
served through the lens, 



/'• 



-\o). 



(1) 



n(6>; > S) 
n (S) 

Here 9 is the position on the sky. In practice, one observes 
the local galaxy distribution down to a limiting magnitude, 
7Ti, so that the cumulative normalized number density is of- 
ten written as 



= M 



L W, 



(2) 



77,(0; < 777) 
no (m) 

with a = dlogiV/dm = 0.4/3 being the slope of the 
(unlensed) differential galaxy number counts. By inverting 
Eq. (El) , it is therefore possible to measure the magnification 
field, fi(d), provided that a sample of background sources 
can be i solated, e.g. on the basis of color and magnitude 
(Sec. 3.4), and that the parameters a and no for this sample 
can be determined to a sufficient accuracy. In the following 
sections, we shall discuss these issues, and present detailed 
analysis of the observed density field n/no in the form of 
radially averaged depletion curves and maximum likelihood 
analysis. 



3.1 Source extractions 

The first step in the analysis of the final images is to ex- 
tract the sources. Source extraction was done with both 



SExtractor (Bertin & Arnouts (1996)) and ImCat (Kaiser 
et al. (1995)), with comparable results. The ImCat results 



were used to identify stars, based on half-light radii, but 
the rest of the analysis below is based on the information 
extracted by SExtractor. The full set of SExtractor object 
catalogs used here is available at the CDS and ADC catalog 
archives. 

We have performed extensive analysis of simulated 
data, using noise characteristics from the final reduced im- 
ages, to find optimal adjustm ents of the va rious parameters 
of the detection software (see Greve (200C) for details). For 



SExtractor, these considerations resulted in choosing a de- 
tection threshold of la per pixel, combined with a Gaussian 
filtering kernel (FWHM of 5 pixels, roughly corresponding 
to our seeing of 6 pixels) , and a minimum detection area of 
10 pixels. The complete SExtractor parameter files used are 
available at the aforementioned website. 

When selecting objects for further analysis from the cat- 
alogs produced by SExtractor, we leave out objects which 
have internal flags larger than 3. This excludes objects with 
saturated pixels, objects that have been truncated, and ob- 
jects with corrupt isophotal or aperture data. 

3.2 Galaxy number counts 

The original motivation for using U band observations for 
depletion studies was the indication of a break in t he galaxy 
number cou nts at faint U magnitudes reported by Williams 
et al. (1996 ). They found a logarithmic slope a « 0.4 for 
£^300 < 26 dropping to 0.05 in the range 26 < C/300 < 28. 
A more detailed analysis of the HDF number counts by 
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Pozzetti et al. (1998J ) shows that the slope drops to a — 0.135 
at C/300 = 25.75, with an allowed range of 0.08 < a < 0.16. 
The differential number counts in our field are shown 
in Fig. hi for the R and U band. The overall higher level 
of br ight R band object s in our data relative to the c ounts 



from Smail et al. (1995 ), Williams et al. (1996), and |Hogg 
et al. ( }997[ ) is simply due to the presence of the cluster. A 
similar effect is seen in the U band. The net result of this 
and the incompleteness at faint magnitudes is to lower the 
logarithmic slope of the differential number counts. This is 
evident from Fig. |lp, where we find a slope of 0.21 in the 
magnitude range 20 <-Rab< 25, compared to an estimate of 
0.27 for the other data in the same magnitude range. 

It is evident that our data is incomplete beyond Uab~ 
26 and 7?ab~ 25.5. This is c onsistent with our completeness 
estimates, but as shown by Gray et al. (2000), the incom- 
pleteness does not affect the depletion analysis, as long as 
the field used to estimate the average background density 
has the same completeness characteristics. It means, how- 
ever, that we are unable to determine the slope at the faint 
end of the luminosity function, and must rely on estimates 
from deeper observations. The slope in the U ma gnitude in- 
terval used in the depletion analysis (section 3.6) is likely to 
be varying fro m a value close to 0.4 at Uab = 24 to the value 
0.135 found by [Pozzetti et al. (1998J ) at Uab~ 26. Similarly, 
the R band slope is dropp ing from a value o f q_r — 0.3 or 



so at bright magnitudes (smail et al. (1995); Hogg et a^ 
(19970 to 0.18 at the faint end (i? A B=27) of the selected 
magnitude interval, in accordance w ith the slopes found i n 
the V a nd I bands from HDF da ta ( |Williams et al. (1996[ ), 
see alsolMadau & Pozzetti (2000J)). 



3.3 Color-magnitude diagrams 

We find the (U — R)ab color by first detecting and measur- 
ing the magnitudes in the U band. The U catalog is then 
matched to the R band image by running SExtractor in "as- 
soc" mode with a search radius of 6 pixels, appropriate to 
compensate for our seeing and error in the alignment be- 
tween the U and R images. In this way, 1229 objects are 
found when associating the U image with the full R catalog. 
These objects comprise our "matched" catalog. 

In Figs. Ha and 0b we show the color-magnitude dia- 
gram for all the 1169 objects in the matched catalog after 
removal of 60 stars (the distinction between stars and galax- 
ies was based on the analysis of half-light radii). The formal 
detection limits in U and R are indicated by dashed lines, 
and the color and magnitude cuts used to select background 
objects are shown with dash-dotted lines (see Sec. 3.4 ). The 



bluing trend at faint R and U magnitudes (e.g. Hogg et al 



(1997)) is reflected in the negative slope of the lower edge of 
the "wedge" shaped region that the objects populate. This 
effect is especially pronounced in Fig. Ha. 

The color- magnitude diagrams clearly show a lower den- 
sity of objects around (U — 7?)ab=3, also seen in U and R 
band data we have obtained for the cluster MS1621+2640 
(van Kampen et al., in preparation). By combi ning our pho- 
tomet ric d ata with redshifts from the survey of |Czoske et al 



distinction is not seen for the field galaxies in the redshift 
survey, suggesting that mostly cluster galaxies will be dis- 
carded by selecting only galaxies bluer than (U — R)ab=3. 
This can also be inferred from the distribution of confirmed 
cluster galaxies, marked by '+', in the color-magnitude dia- 
grams. The radial gradient observed fo r g — r color in many 
clusters (see e.g. Morris et al. (1998)) is also seen in the 
(U — R)ab color. These features are suggestive of a strong 
correlation between star formation activity in cluster galax- 
ies and the cluster environment. We shall investigate this 
aspect of the data in greater detail in a forthcoming paper 
(van Kampen et al., in preparation). 



3.4 Selection of background candidates 

For the depletion analysis, it is important to remove cluster 
and foreground galaxies from the object catalogs as the num- 
ber density of cluster galaxies typically increases towards the 
center, interfering with radial gradients in the number den- 
sity of the background. In line with the discussion above on 
color-magnitude diagrams, we chose to use simple limits on 
magnitudes and color, restricting the depletion analysis to 
the magnitude ranges 24 < (R, U)ab < 27 and to objects 
bluer than (U — R)ab = 3. The limiting U band magnitude 
of 27.2 implies that R band objects in the magnitude range 
24 < Rab < 27 that have a well defined counterpart in 
the U band image are all bluer than (U — -R)ab=3, as can 
be seen from figure Ha. Therefore, the exclusion of objects 
redder than (U — -R)ab=3 will only affect the selection of 
background objects in the U band. 

The background samples are constructed by first se- 
lecting objects in the magnitude ranges 24 < (R, U)ab < 27 
from the full R and U band catalogs. In order to exclude 
the red (cluster) objects, one could in principle extract the 
background sample from the matched catalog, with the ad- 
ditional constraint on color. Doing this, however, one tends 
to miss primarily the faintest objects, partly because some 
of the U objects are simply too faint in R to be detected, 
and vice versa, as well as because the association is more 
difficult the fainter the object is. We avoid this drawback by 
constructing masks (using the shape parameters estimated 
by SExtractor) covering the (U — R)ab > 3 objects in the 
matched catalog. The objects that fall on these masks are 
then deleted from the background samples. 

In addition, the masks are augmented with the bright 
objects ((R, U)ab < 24), to cover the parts obscured by 
the foreground and cluster candidates in the correspond- 
ing image (see figure W). The estimated (Poisson) errors for 
the depletion curves are corrected for the effective area of 
each annulus of the radial averaging by subtracting off the 
number of pixels in the mask from the pixelized area of the 
annulus. The correction due to the mask is quite important 
for the innermost bins, especially in the R band, where the 
bright cluster galaxies obscure a large fraction of the area. 
Hence, the innermost bins have the highest uncertainty due 
to the small effective area. 



(1999) and Dressier et al. (1999) for the CL0024+1654 field, „ _ „ , ,. . . 

— ca u il — 777 — m — — i — c u **•" Depletion analysis 



we find that the (U — R)ab color of the cluster galaxies shows 
a distinct bi-modal distribution, with a minimum around 
(U — R)ab=3 and a large fraction above (U — R)ab=3. This 



Radial averaging is typically used to boost the significance 
of the depletion signal (e.g. Taylor et al. (1998); Mayen & 
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Figure 1. (a - ) The Ra w number counts compared to i? band data from Hogg et al. (1997), and Smail et al. (1995), and I band data from 
Williams et al. (1996). The effect of the cluster is mor e pronounced at brigh t mag nitudes, while inco mpleteness sets in around i?AB=25. 



(b) The [Tab number counts compared to data from Williams et al. (199C), and Hogg et al. (1997). Again, the presence of the cluster 
results in overall higher counts, and the turn-over at faint magnitudes reflects the incompleteness beyond Uab = 26. 
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Figure 2. (a) The (U — R)ab vs Rab color-magnitude diagram. The formal detecion 3<r limits, t/AB,lim= 27.2 mag and RAB,lim= 27.5 
mag, are indicated with dashed lines. Dash-dotted lines indicate the color (([/ — -R)ab< 3) and magnitude (i?AB> 24) limits used to 
mask out foreground objects in the R band. The upper edge of the "wedge" shape is defined by the limiting magnitude in U, while the 
slope of the lower edge clearly reflects a bluing trend towards fainter magnitudes. Note also the distinct low density around (U — -R)ab=3. 
Confirmed cluster galaxies are marked by '+'. (b) The (U — R)ab vs t^AB color-magnitude diagram. Dash-dotted lines show the color 
((U — -R)ab5= 3) and magnitude (Uab> 24) limits used to mask out foreground objects in the U band. The general features are the same 
as in (a) 



Soucail (2000)). The number density is estimated by count- 
ing the number of objects in annuli around the cluster cen- 
ter. This average density is normalized by an (indepen- 
dent) estimate of the unlensed background density of the 
same class of objects resulting in a radial depletion curve. 
The width of the annuli, Ar, may be adjusted to obtain a 
reasonable balance between smoothing and shot noise. We 
calculate the number density averaged over 10" wide an- 
nuli as a function of radius from the cluster center, taken 
to be the center of light (a(J2000) = 00 h 26 m 35.8 s and 
S(J2000) = +17°09'42.3") of the R image within 100" ra- 
dius after a few bright stars and a foreground galaxy have 
been masked out. The number counts are normalized to 
the background density, no, estimated as the average den- 
sity of the corresponding objects in the annulus between 



150" < r < 240". This assumes that the cluster does not 
affect the background significantly beyond 150". (Ideally, 
the unlensed background density should be estimated from 
a separate field, known not to be affected by strong fore- 
ground lenses.) We find values of no(U) — 28 ± 2 arcmin -2 
and fio(-R) = 49 ± 7 arcmin -2 , with 1084 and 1780 selected 
background objects in the U and R band, respectively. (The 
uncertainties quoted are jackknife estimates, see below.) 

The obscuration by the masks (separate for each band) 
is taken into account in determining an effective area, 
A c s(ri), for each annulus, n — Ar/2 < r < n + Ar/2. The 
uncertainties may as a first approximation be taken to be 
due only to Poisson noise in the counts, so that 

AN(n) = 7Vcx P (r0 1/2 = {n A cB (r t )) 1/2 
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Figure 3. (a) The R band mask used in the depletion analysis superimposed on the R band image of CL0024+1654. Objects brighter 
than -Rab= 24 and redder than (U — R)ab = 3 are masked out. The background objects are marked with '+'. (b) Same as (a) for the U 
band. Objects brighter than Uab= 24 and redder than (U — R)ab= 3 are masked out. 



where N(ri) are the actual number counts in the annulus 
centered on n, and -/V oxp (ri) — noA c s(ri) is the expected 
number of counts in the annulus. Then, 



An(n 
n 



(n A eff (r;)) 



-1/2 



Given an uncertainty in the absolute background density, 
Ano, a more conservative estimate of the uncertainty is 



n(rj) 
no 



An(r.) 



no 



+ 



n(rj)Ano V 
nl ) 



1/2 



(3) 



We have made simple estimates of the intrinsic uncertainty 
in no due to non-uniformity of the background, effects of 
bright objects, etc, by masking out each of the four quad- 
rants of the images in turn, numbered counter-clockwise 
from the NE quadrant. The four values of the background 
obtained in this way are then used to form a standard jack- 



in Fig. ki A clear sign of depletion is seen in both bands. In 
the R band the density falls steadily from a radius of around 
100" all the way to the center. Most lens models predict 
a turn-up in the density inside the innermost critical line; 
the indication seen here in the innermost two bins is non- 
significant due to the obscuration of the foreground galaxies. 
On the other hand, it is easier to detect faint background 
objects closer to the cluster center in the U band, since the 
dominant cluster galaxies are all quite red. Given that we 
are successful in identifying cluster members, the observed 
turn-up is quite sharp and significant in the U band. 

Our data is not complete to sufficiently faint magni- 
tudes to allow estimation of the slopes an and air, for the 
magnitude ranges involved. However, we find that most of 
the depletion signal in the U band is due to galaxies with 
24 < Uab < 26.5, which indicates that the slope has already 
dropped well below 0.4 i n this magnitude ra nge. This is at 
odds with the results of Hogg et al. (1997), who found a 



slope of 0.467 down to U13 = 25.5, roughl y corresponding 



knife estimate the variance of no. We find that or/ (no) = 2.2 

and n a [n. n ) = 7.2. The errors indicated on Fig. §a andgh in- to C/ab=26.3. The recent investigation by [Volonteri et al 



elude a lfractional uncertainty n = Ano /no of nn — 7.2/49 = 
0.15 and rju = 2.2/28 = 0.08, respectively. Note, that an 
error in no only affects the normalization of the depletion 
curve, while the relative depletion does not change. The es- 
timated background density in the R band is strongly influ- 
enced by quadrant 2, where the density of selected objects 
is markedly higher than in the other three quadrants. The 
form of the depletion curves is however robust. The varia- 



(20001) of the HDF-S shows no evidence of a break either, 



down to the faintest magnitude bin of C/ab=26.5. 

Although the data may not warrant a detailed quanti- 
tative interpretation of the depletion curves, especially when 
taking the uncertainty in the slope, a, and unlensed back- 
ground estimate, no, into account, it is in principle possi- 
ble to use multiband depletion measurements to constrain 



tion in the R band may at least partly be explained by a 



the properties of the background population (se e e.g. Tay- 



higher aackground flux in quadrants 3 and 4, due to bright 



stars. The background flux in the U band is more uniform 
across the whole field. The estimated values of n should only 
be taken as an indication of the level of uncertainty in no. 
The depletion curves in the R and U band are shown 



lor et al. (1998J ); |Bartelmann fe Schneider (1999J ); |Mayen & 
Soucail (2000p , and references therein) due to the depen- 
dence on the slope of the luminosity function. The redshift 
distribution of the background objects will also affect the 
width and position of the minimum of the depletion curve, 
essentially through the dependence of the Einstein radius on 
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the lens-source distance. In Fig. H, we compare the radially 
averaged data to the sing ular isothermal sphere (SIS) model 



described in section 3.6. In case of the R band depletion 
curve, results for an — 0.15 and 0.2 are shown, with the 
Einstein radius te = 25". In addition, we show the theo- 
retical depletion curve for te = 40" and an — 0.32. The 
U band depletion is compared to models with ajj — 0.15 
and 0.2, with a fixed Einstein radius of te = 20", as well 
as to ctu = 0.32 and te = 32". (This choice of parameters 
is dictated by results of the maximum likelihood analysis 
in section 3.6.) We note that the depletion curve is better 
represented by the te = 25" and te 
and U bands, respectively. 



20" models in the R 



3.6 Einstein radius estimates 

We use the maximum lik elihood approach developed by 
Schneider et al. (2000 ) and pray et al. (2000 ) to investigate 
the depletion in more detail. A singular isothermal sphere 
(SIS) model, where the magnification is a simple function of 
the distance from the center of the lens, r, and its Einstein 
radius, 



Msis(r;r E ) = 



-r E /r 



is used to calculate the likelihood as a function of a single pa- 
rameter, the Einstein radius te- Recall, that the theoretical 
depletion curve is in this case given by 



n(r) 
n 



1 



1 - r E /r 



where the slope of the luminosity function is /3 = 2.5a. For 
simplicity, we adopt a constant slope, noting that a proper 
treatment of our data would take the variation in a into 
account, since we are very likely probing a transition region 
in both the R and U bands. We take the (unknown) uncer- 
tainty in the unlensed backgro und density into a ccount, us- 
ing the likelihood function (see Gray et al. (2000), Eq. (10)), 



log L™=lP = -noI+{0-l)J2lDgtiOi 



(n - n ) 2 
2(nn ) 2 



(4) 



(•>) 



+iV log no - 

where I = J d 2 0(u(0))' 3 ~ 1 , and 

— = v + [y + n N) ' 
no 

with v — ^(1 — nnol). The first term in Eq. (kf) stems from 
the likelihood of observing a total of N galaxies in the whole 
field, while the second term is the log-likelihood for finding 
those N galaxies at the observed positions 6i = (ri,8i) in 
polar coordinates. The third term is just a normalization, 
but is needed in this case as we include prior information 
on no, taken to be a Gaussian with mean no and dispersion 
nno. The maximum likelihood estimate of no is then given 
by Eq. (|5j), and is inserted back into Eq. (||) to obtain the 
log-likelihood. Note, that the model parameters (te in the 
SIS case) are all hidden in the magnification u. 

To assess the effect of the uncertainty in no we study 
four values of the dispersion by setting n = 0, 0.05, 0.1, and 



0.15. The highest value is probably a conservative estimate 
taking clustering of the background into account, given the 
l atest measurements of the two point correlation function 
(Fynbo et al. (200C)), which indicate insignificant correla- 
tion beyond 10 for the magnitude range considered here. 
The estimated amplitude of the tw o point correlation fun c- 
tion in the U band is even lower (Brunner et al. (200C)), 



making the U band depletion less sensitive to errors due to 
clustering. Possible bias in the estimated unlensed density 
is studied by shifting no up by 10%. We would expect the 
estimated no to be biased low relative to the true value no, 
since the depletion due to the cluster probably extends be- 
yond the inner radius n = 150". Contamination by cluster 
members would however work in the other direction, and it 
is difficult to disentangle the two effects with our data. By 
varying the inner radius n from 150" to 220" we find that 
the U band background density is quite stable, while in the 
R band no does indeed increase by roughly 10%, going to 
54 arcmin" 2 at n = 200". 

We also study the dependence of the estimated Ein- 
stein radius on the slope in the range 0.1 < a < 0.4. In the 
U band, au = 0.135 is appropriate for Uab> 25.5, but since 
50% of the selected background galaxies in the magnitude 
range 24 <(7ab< 27 are brighter than this, we expect the 
effective slope to be steeper. Similarly, in the R band, the 
slope at the brig ht end is of order 0. 32, dropping to ~ 0.1 8 



at the faint end (Small ct al. (1995 ); Williams et al. (1996 ); 



Madau fc Pozzetti (2OO0| )). As 50% of the background can- 



didates have Rab> 25.3, the effective slope probably lies in 
this range. The uncertainties, ArE, quoted below are esti- 
mated as the width of the likelihood function, LjJ 1 , 



ir E 



J(r E -r E fL^dr^ 
fL^dr> E 



(6) 



where te is the estimated value of the Einstein radius, cor- 
responding to the maximum of L^ . 

The results of this analysis for the R and U band are 
shown in Fig. 0, where the estimated Einstein radius, te, is 
plotted as a function of the slope a for values of the disper- 
sion parameter n in the expected range — 0.15. The error 
bars indicate ±ArE- The estimated te grows with increasing 
a, up to a — 0.4 where the theoretical depletion signal van- 
ishes. The effect of increasing a is to lower the log-likelihood 
at small te, so that local maxima at higher te take over as 
global maxima, causing discrete jumps in the (a, te) plots. 
Increasing n tilts the I J 1 curve upwards at high te, but for 
n > 0.05 the estimated te is not very sensitive to this param- 
eter. Assuming that our unlensed background estimate no is 
biased low relative to the true value by about 10%, we find 
that even the n = estimates of te are compatible with the 
n > 0.05 cases. The dependence of the log-likelihood IjJ 1 on 
the parameters a and n is shown in Fig o where we plot Ij 1 
as a function of the Einstein radius, te, for selected values of 
these parameters. We simply tabulate l}] 1 ' with ArE = 0.5" 
and take the estimated Einstein radius to be the te where 
the tabulated function is maximized. The dips in the Ij 1 
curves are caused by galaxies that happen to lie close to the 
current te (the critical line being tested) , where the second 
term in Eq. H dominates, and the probability of finding a 
galaxy vanishes. 
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Figure 4. (a) Normalized density of faint R band galaxies as a function of radius from the center of CL0024+1654. The selected 
magnitude range is 24 < -Rab < 27, and the radial binsize is 10". The upper radial scale shows the distance from the center in h^ 
kpc, assuming an Q = 1 cosmology. The data is compared to three singular isothermal sphere (SIS) models, with {rE,ct}={25" , 0.15}, 
{25",0.2}, and {40", 0.32}. (b) Same figure for the U band. The selected magnitude range is 24 < U A B < 27. Here, the SIS models shown 
have {r E ,a}={20", 0.15}, {20" ,0.2}, and {32" ,0.32}. 
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Figure 5. (a) The estimated Einstein radius in the SIS model as a function of assumed slope a for the four values of the dispersion 



parameter r\ = 0, 0.05, 0.1, and 0.15, using the R band depiction. The biased model has no 
kpc, assuming an Q = 1 cosmology, (b) Same as (a) using the U band depletion. 



l.lno and n = 0. The right scale is in h 5 



From Fig fja we see that for slopes a < 0.25 the es- 
timated Einstein radius is 25" ± 3" in the R band with a 
limited range of parameters suggesting a lower value of 20" . 
The unbiased r\ = case results in yet a lower value of 
around 13". For the rest of the a range, the maximum of 
Ijl 1 occurs around 40", again with the unbiased r\ = case 
as an exception. Similar features are seen in the U band 
(Fig. Bp) , except that for any given set of parameters the te 
estimates are always lower than in the R band. In any case, 
the U band Einstein radius of 20" ± 3" found for a < 0.25 
and 0.1 < r\ < 0.15 is within the uncertainty of the other 
estimate of 17" ± 3" found for the r\ — 0.05 and the biased 
models The unbiased r] = fl model is compatible 



'I 



for a < 0.28 where the n?, 



17" peak in W 



with te_ 

takes oyer. Since the giant arcs seen at around 30" distance 



from trie clus ter center are the multip le images of a. galaxy 
at z = 1.675 (Broadhurst et al. (2000)) it seems reasonable 



to reject any estimates that exceed this value. It must how- 
ever be emphasized that this analysis is based on a single 
underlying model (SIS) and the final constraints on the pa- 
rameters a, no, and n must come from reliable deep galaxy 
counts. To summarize, our analysis indicates that the Ein- 
stein radius may be slightly lower in the U band compared 
to the R band, suggesting a lower median redshift for the U 
band objects. 

For comparison, we have also calculated the likelihood 
function for a truncated isothermal sphere (TIS) model, with 
a core radius of r c — 10" in accordance with results found 
from s trong lensing ( [Tyson et al. (1998|); Broadhurst et al 
(2000| )) and X-ray observations ( |Bohringer et al. (200^TT 



The T IS model is given by equation (8) in Mayen fc Soucail 



(2000| ) 



The estimated Einstein radius in the TIS model is 
almost identical to the SIS results, as expected from the 
small core radius. Finally, we have checked that our re- 
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Figure 6. (a) Maximum likelihood estimation of the Einstein radius in the SIS model using the R band depletion. The log-likelihood is 
plotted for the same parameters as the theoretical depletion curves in Fig. kja, i.e. r? = 0.1 and a = 0.15, 0.2, and 0.32. In addition, the 
log-likelihood for the TIS model (see text) with r\ = 0.1 and a = 0.15 is shown. The upper radial scale shows the Einstein radius in h^ 
kpc, assuming an Q = 1 cosmology, (b) Same as (a) using the U band depletion. 



suits ar 3 insensitive to a more stringent selection of R band 



objects by redoing the analysis for the magnitude range 
25.3 <-Rab< 27. For a < 0.2 and rj > 0.05 this shifts the 
estimated Einstein radius slightly outwards, to around 27". 
This is to be expected, as the fainter sample will have a 
higher median redshift. 



4 DISCUSSION AND CONCLUSIONS 

The applicability of lensing methods to determine cluster 
masses and density profiles has been a major issue since 
it became viable to detect lensing signatures. Using strong 
lensing effects, where background galaxies are stretched into 
giant arcs and/or show up as multiple images, it is possible 
to probe the cluster potential roughly out to the Einstein ra- 
dius. On larger scales, the shear and depletion due to the lens 
may be used to constrain its mass distribution, once these 
lensing effects are disentangled from the intrinsic properties 
of the sources. 

From simple arguments, it can be shown that shear 
detections have a higher signal-to-noise ratio than magni- 
ficatio n detections with the number cou nt method (Mellier 
(1999| ); |Bartelmann fc Schneider (199SJ )) 



(2000D ). 



It must however 
be emphasized, that from the observational point of view 
it is easier to obtain deep galaxy number counts than de- 
tailed shear maps that require good seeing and accurate de- 
termination of the PSF for deconvolution. Note also, that 
ground based shear detections are limited to galaxy sizes 
above the seeing limit. With 0.5" seeing in the R band this 
corresponds to Rab< 25, giving a galaxy number density 
of roughly 30 arcmin -2 . Depletion analysis on the other 
hand benefits from detections up to the magnitude limit, so 
that with deep enough observations one can get much higher 
counts. Another aspect in favor of the number count method 
is that it is not subject to the "sheet-mass d egeneracy" of 
the shear analysis (e.g. Schneider et al. (200C )). This does, 
however, not come for free, since absolute surface density 
measurements using depletion require an a ccurate determi- 
nation of the unlensed background density (pchneider et al 



The intrinsic uncertainties due to the combination 
of observational errors, including contamination from fore- 
ground objects and false detections, and clustering of the 
background population, certainly plague the number count 
method. In some wavelength bands these uncertainties make 
it diffi cult to obtain sign ificant constraints on model param- 
eters ( |Gray et al. (200C| )). 

Our data and analysis show that detailed inferences 
from the present data are to some extent sensitive to poorly 
known a priori parameters (a, no, r), ...) making it dif- 
ficult to set strong limits on the parameters of the lens 
model. However, we do find compatible results for the two 
bands when the slope of the differential number counts of 
the background objects, cv, is lower than 0.25 independent of 
the uncertainty in the unlensed background densities. It is 
worth noting, that the U band depletion curve is both shal- 
lower and nar rower than the R band d epletion, suggesting a 
steeper slope (Mayen & Soucail (200C)) in the U band. The 
two extreme values of the Einstein radius, r^ = 17" and 
te = 25", found from the U and R band data respectively 
for a < 0.25, may be used together with the observed Ein- 
stein radius, rE, ga i, a ssociated with the singl e background 
galaxy at z = 1.675 ( Broadhurst et al. (2000 )) to estimate 
the median, (zs), of the redshift distribution of the back- 
ground galaxies. In the case of a flat £1 = 1 cosmology, we 
find that the ratio rE/VE, ga i = 17/30 ~ 0.57 corresponds 
to {zs) = 0.7 while 7"E/»"E,gai = 25/30 ~ 0.83 corresponds 
to (zs) = 1.1. This latter value is in good agreement with 
the models of Mayen & Soucail (200C), although the incom- 



pleteness of our data at faint magnitudes and the different 
filters used make a direct comparison difficult. Note also, 
that taking into account the uncertainty on the typical U 
band Einstein radius estimate of te{U) = 20" ± 3" and the 
R band estimate of te(-R) = 25" ± 3" these two values do 
not differ significantly. 

The maximum likelihood analysis estimates of the Ein- 
stein radius using a singular isothermal sphere model is con- 
sistent with the observed giant arcs, provided that the slope 
a < 0.25. The estimated Einstein radius is independent of 
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the dispersion of the unlensed background density, parame- 
terized by r/, provided that our estimate of the background 
density is biased low. We would indeed expect that to be 
the case since the cluster is likely to deplete the background 
density beyond the 150" assumed in our estimation of the 
background density. We conclude that although it is difficult 
to draw strong conclusions about e.g. the different galaxy 
populations probed with the two bands without a more rig- 
orous treatment of the effects of a variable a and firmer 
estimates of the unlensed background densities, we do get 
results that are consistent with results obtained by different 
methods by others. 

One of the main aims of our work was to look for de- 
pletion of background galaxies in the U band. Our data 
clearly show this effect towards the cluster CL0024+1654, 
giving independent evidence for a flattening to a slope /3 < 1 
(q < 0.4) at the faint end of the U band luminosity func- 
tion. While the reality of this flattening is still debatable 



(e.g. Pozzetti et al. (1998); Volonteri et al. (2000)), further 



depletion studies in the U band will help to constrain the 
slope of the number counts in this band, ft is clear that ul- 
timately the shape of the number counts at faint U band 
magnitudes must be determined from deep, reliable obser- 
vations in the field. This would also settle the absolute un- 
lensed density of background sources, both of which is nec- 
essary for making U band depletion a reliable tool in cluster 
studies. 

As the dominant cluster galaxies are relatively faint in 
the U band, we are able to measure the background number 
density close to the cluster center accurately enough to find 
a significant turn-up in the density as expected from simple 
lens models. This strengthens the indication of a turn-up 
in the R band, where the estimated background density is 
starting to increase before reaching the innermost parts of 
the cluster where the foreground objects cover a too large 
fraction of the available area to warrant a significant mea- 
surement. 

With the present data, where the background density 
is of order 30 — 50 arcmin -2 , the inherent uncertainties of 
the number count method still give rather loose bounds on 
lens model parameters, both when fitting models to radially 
averaged depletion curves, and when more elaborate maxi- 
mum likelihood analysis is applied. It goes without saying 
that with deeper number counts towards lenses, and bet- 
ter estimates of the properties of the unlensed background 
population, the reliability of the number count method will 
increase. Some of the difficulties of the parametric modeling 
inherent to these methods may however also be alleviated 
by using more sophisticated, non-parametric methods. We 
plan to study this in more detail for the data presented in 
this paper. 



ACKNOWLEDGEMENTS 

This work was supported in part by the Danish Research 
Foundation (through its establishment of the Theoretical 
Astrophysics Center), The Icelandic Research Council, The 
Research Fund of the University of Iceland, The Danish Nat- 
ural Science Research Council (SNF), and by the "Access to 
Large- Scale Facilities" program of the European commis- 
sion of the EU. The Nordic Optical Telescope is operated 



on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway and Sweden, in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias. We are indebted to Ian Smail, David Hogg, 
and Lucia Pozzetti for providing data on galaxy number 
counts, and to Oliver Czoske for providing redshift data for 
CL0024+1654. We thank the referee for useful comments 
and suggestions. 



REFERENCES 

Athreya R., Mellier Y., van Waerbeke L., Fort B., Pello R., 

Dantel-Fort M., 2000, A&A, submitted 
Bartelmann M., Schneider P., 1999, Weak gravitational 

lensing 
Bertin E., Arnouts S., 1996, A&AS, 117, 393 
Bezecourt J., Kneib J.-P., Soucail C, Ebbels T. M. D., 

1999, A&A, 347, 21 

Bohringer H., Soucail C, Mellier Y., Ikebe Y., Schuecker 
P., 2000, A&A, submitted 

Bonnet H., Mellier Y., Fort B., 1994, ApJ, 427, L83 

Broadhurst T., Huang X., Frye B., Ellis R., 2000, ApJ, 534, 
L15 

Broadhurst T. J., Taylor A. N., Peacock J. A., 1995, ApJ, 
438, 49 

Brunner R. J., Szalay A. S., Connolly A. J., 2000, ApJ, 
accepted 

Clowe D., Luppino G. A., Kaiser N., Henry J. P., Gioia 
I. M., 1998, ApJ, 497, L61 

Czoske O., Soucail G, Kneib J.-P., Bridges T., Cuillandre 
J.-C, Mellier Y., 1999, in Gravitational lensing: Recent 
progress and future goals. A wide-field spectroscopic sur- 
vey in the cluster lens C10024+17. Boston University 

Davis L., , 1999, private communications 

Dressier A., Gunn J. E., 1992, ApJS, 78, 1 

Dressier A., Gunn J. E., Schneider D. P., 1985, ApJ, 294, 
70 

Dressier A., Smail I., Poggianti B. M., Butcher H., Couch 
W. J., Ellis R. S., Oemler A. J., 1999, ApJS, 122, 51 

Dye S., Taylor A., 1998, MNRAS, 300, L23 

Dye S., Taylor A. N., Thommes E. M., Meisenheimer K., 
Wolf C, Peacock J. A., 2000, MNRAS, submitted 

Fort B., Mellier Y., Dantel-Fort M., 1997, A&A, 321, 353 

Fruchter A. S., Hook R. N., 1997, Astron. Data Soft- 
ware and Systems, VI, 147 

Fukugita M., Shimasaku K., Ichikawa T., 1995, PASP, 107, 
945 

Fynbo J. U., Freudling W., Moller P., 2000, A&A, 355, 37 

Gray M. E., Ellis R. S., Refregier A., Bezecourt J., McMa- 
hon R. G, Beckett M. G, Mackay C. D., Hoenig M. D., 

2000, MNRAS, submitted 

Greve T., 2000, Master's thesis, University of Copenhagen 
Hoekstra H., Franx M., Kuijken K., 2000, ApJ, 532, 88 
Hogg D. W., Pahre M. A., McCarthy J. K., Cohen J. G, 
Blandford R., Smail I., Soifer B. T., 1997, MNRAS, 288, 
404 
Hook R. N., Pirzkal N., Fruchter A. S., 1999, As- 
tron. Data Software and Systems, VIII, 337 
Joffre M., Fischer P., Frieman J., McKay T., Mohr J. J., 
Nichol R. C, Johnston D., Sheldon E., Bernstein G, 2000, 
ApJ, 534, L131 



Depletion due to the cluster lens CL0024+1654 11 



Kaiser N., Squires G., Broadhurst T., 1995 

Koo D. C, 1988, in Rubin V. G., Coyne G. V., eds, Large 

Scale Motions in the Universe Princeton University Press 
Landolt A. U., 1983, AJ, 88, 439 
Luppino G. A., Kaiser N., 1997, ApJ, 475, 20 
Madau P., Pozzetti L., 2000, MNRAS, 312, L9 
Mayen C, Soucail G., 2000, A&A, submitted 
Mellier Y., 1999, ARA&A, 37 
Morris S. L., Hutchings J. B., Carlberg R. G., Yee H. K. C, 

Ellingson E., Balogh M. L., Abraham R. G., Smecker- 

Hane T. A., 1998, ApJ, 507, 84 
Mutchlcr A., Fruchter A. S., 1997, Am. Astron. Soc. Meet- 
ing, 191, 4107 
Pozzetti L., Madau P., Zamorani G., Ferguson H. C, 

Bruzual A. G., 1998, MNRAS, 298, 1133 
Schneider D. P., Dressier A., Gunn J. E., 1986, AJ, 92, 523 
Schneider P., King L., Erben T., 2000, A&A, 353, 41 
Shapiro P., Iliev I., 2000, ApJ, submitted 
Smail I., Dickinson M., 1995, ApJ, 455, L99 
Smail I., Hogg D., Yan L., Cohen J., 1995, ApJ, 449, L105 
Soucail G., Ota N., Bohringer EL, Czoske O., Hattori M., 

Mellier Y., 2000, A&A, 355, 433 
Taylor A. N., Dye S., Broadhurst T. J., Benitez N., van 

Kampen E., 1998, ApJ, 501, 539 
Tyson J. A., Kochanski G. P., dell' Antonio I. P., 1998, ApJ, 

498, L107 
van Kampen E., 1998, MNRAS, 301, 389 
Volonteri M., Saracco P., Chincarini G., Bolzonella M., 

2000, A&A, accepted 
Williams R. E., Blacker B., Dickinson M., Dixon W. V. D., 

Ferguson H. C, Fruchter A. S., Giavalisco M., Gilliland 

R. L., Heyer I., Katsanis R., Levay Z., Lucas R. A., MCE1- 

roy D. B., Petro L., Postman M., Adorf H. M., Hook R., 

1996, AJ, 112, 1335 



